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We review recent one and two loop MS Landau gauge calculations using the Gribov-Zwanziger Lagrangian. 
The behaviour of the gluon and Faddeev-Popov ghost propagators as well as the renormalization group invariant 
effective coupling constant is examined in the infrared limit. 
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I. INTRODUCTION. 

Quantum Chromodynamics (QCD) is the quantum field 
theory describing the physics of the strong nuclear force. In 
the ultraviolet regime the theory describes asymptotically free 
quarks and gluons, |l|, |2y, which at high energies behave as 
fundamental particles with propagators of the form 1 / p 2 as 
the momentum p increases. However, free quarks and glu- 
ons are not observed in nature as isolated units due to the 
confinement mechanism. Moreover, there is a wealth of ev- 
idence to indicate that at low energies the behaviour of the 
gluon propagator differs significantly from that in the ultra- 
violet. (For example, see the recent reviews of y| |4j].) For 
instance, both Dyson Schwinger equations (DSE) and lattice 
regularization observe that the gluon propagator vanishes at 
zero momentum, yl la LZl l2l 1 1 Qfl . This property is referred 
to as gluon suppression. Whilst it has been established in the 
Landau gauge several other features emerge in the infrared. 
The propagator of the associated Faddeev-Popov ghost field 
diverges at zero momentum more rapidly than its l/p 2 ul- 
traviolet behaviour. Again this ghost enhancement has been 
observed in DSE and lattice studies where both approaches 
obtain good agreement on the value of the exponent in the 
power law form of the propagator, JH 0, 0, H] • Another in- 
frared property which is not fully resolved by DSE and lattice 
analyses is that of the freezing of a particular definition of the 
effective strong coupling constant. Whilst it is accepted that 
the particular quantity does not diverge at zero momentum, it 
is not clear whether it freezes to a finite or zero value. Given 
that these infrared properties of Landau gauge QCD have been 
well established by the numerically intensive DSE and lattice 
regularization techniques, one natural question to pose is, can 
such phenomena be probed from a Lagrangian point of view, 
and if not initially quantitatively, can it be accessed and under- 
stood qualitatively? This is therefore the purpose of this arti- 
cle, which is primarily to review recent activity in this area, 
iflll [131 . by using Zwanziger's elegant study on the Gribov 
problem, Eld Eld- 

Briefly, Gribov pointed out in [ 17] that one cannot uniquely 
fix the gauge globally in a non-abelian gauge theory. Specif- 
ically for some gauge configurations one can construct one 
or more copies which can be related by a gauge transfor- 
mation. To resolve this Gribov proposed that the path inte- 
gral be restricted to the region defined by the first zero of the 
Faddeev-Popov operator which is defined as the Gribov vol- 
ume. Whilst one can still have gauge copies inside the Gribov 



horizon, CI EllIS, within it there is a smaller region known 
as the fundamental modular region where there are no Gribov 
copies. The key feature of the restriction to the first Gribov re- 
gion is the introduction of the Gribov mass scale, y, which can 
be related to the volume of the first Gribov region. The restric- 
tion of the path integral significantly modifies the behaviour 
of the gluon and ghost propagators in the infrared, I17I1 . For 
instance, the (non-local) Gribov Lagrangian contains a gluon 
which is infrared suppressed and an enhanced Faddeev-Popov 
ghost. The latter derives from the one loop ghost 2-point func- 
tion where the gap equation satisfied by the Gribov mass is 
central to the one loop enhancement. Whilst Gribov's Lan- 
dau gauge analysis gave significant insight into the infrared 
properties, one difficulty was the inability to easily perform 
practical computations since Gribov's approach effectively 
amounted to using a Lagrangian with a non-locality. The 
breakthrough in this area came with Zwanziger's construction 
of a localised Lagrangian in |13lll4l 15 LI16I1 . Significantly that 
Lagrangian is renormalizable, ll2ll l22ll . and therefore allows 
one to carry out calculations. 



II. GRIBOV-ZWANZIGER LAGRANGIAN. 

The Gribov-Zwanziger Lagrangian involves two pieces, 
0]. One is the usual Landau gauge fixed QCD Lagrangian 
and the other part involves additional fields which we will 
refer to as Zwanziger ghosts. These are {§^ b ,§^ b ,03 c ^ h ,03^'} 
where the latter pair are anticommuting. The full localised 
renormalizable Lagrangian is, (u 



l qcd + pbnp ( Dv ^y<b _ ffbiqi (DvC0/j ) 
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dN A f 
2g 2 



(1) 



where L@ CD denotes the usual Lagrangian involving the gluon 
and Faddeev-Popov ghost fields c" and c", is the usual 
covariant derivative, f abc are the colour group structure func- 
tions, Na is the dimension of the adjoint representation, g is 
the coupling constant, and d is the spacetime dimension. The 
commuting Zwanziger ghosts and implement the Gri- 
bov horizon condition which in the Gribov formulation is de- 
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fined by the vacuum expectation value, [17], 
The equivalent definition in ([TJ is, (3, 



(2) 



(3) 



and these two equations define the gap equation satisfied by 
the Gribov mass y. The relation between both definitions of 
the horizon conditions is made manifest from the equation of 
motion of the Zwanziger ghost 



cctbc 



(4) 



Further, focusing on the part of the Lagrangian involving 
purely gluons and using then the non-locality of the origi- 
nal formulation of Gribov becomes apparent 



L = — X -G%G a ^ 
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1 
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(5) 



As with any Lagrangian one needs to determine the explicit 
values of the renormalization constants and hence anomalous 
dimensions in order to perform practical calculations. As was 
demonstrated in Il2ll 12211 the Landau gauge Gribov-Zwanziger 
Lagrangian has an interesting renormalization structure. In 
addition to the renormalization constants of the usual QCD 
Lagrangian the new fields and parameters each have a renor- 
malization constant 



4>o 



ab/i 



CO, 



where the subscript o denotes the bare quantity. However, it 
was demonstrated in 12 ll 12211 by the algebraic renormaliza- 
tion programme that to all orders in perturbation theory in the 
Landau gauge 



and 



Z c — Za — z a 



7 _ 7 l/4 7 l/4 

ZVy — ZV^ ZV C 



(7) 



(8) 



where Z A and Z c are respectively the gluon and Faddeev- 
Popov ghost renormalization constants. Hence, no new renor- 
malization constants are required and, moreover, the original 
QCD anomalous dimensions for the gluon, Faddeev-Popov 
ghost and quark as well as the p-function remain unchanged. 
So the presence of the localising fields does not alter the ul- 
traviolet running. For completeness we note that the anoma- 
lous dimension of y is known to three loops in MS for arbi- 
trary colour group and at four loops for SU (N c ) in the Landau 
gauge, iH. 

Whilst this article primarily deals with the Landau gauge, 
as an aside one can begin to address what happens in other lin- 
ear covariant gauges aside from the Landau gauge by restor- 
ing the usual a dependent part to L@ CD in ((TJ. In I2411 the 



construction of the horizon function has been examined for 
small a where a is the gauge fixing parameter in a linear co- 
variant gauge. In this approximation one can repeat some of 
the Landau gauge construction and establish a gap equation. 
In this context one can examine the renormalization of (Q]) 
where 1_Q CD is replaced by the a dependent Lagrangian. It is a 
straightforward exercise to examine the three loop renormal- 
ization for non-zero a which has not been determined before. 
The machinery to do this involves the application of the MIN- 
CER algorithm, [25, 26], written in the symbolic manipulation 
language Form, 12711 . and the intensive use of computer al- 
gebra. The full three loop renormalization of the a extended 
Gribov-Zwanziger Lagrangian has been performed in the MS 
scheme and the following a dependent anomalous dimensions 
emerged for the additional fields and y to three loops 



lc{a) 
J y (a) 



y<j,(a) = y w (a) 
[16T F Nf-{35 + 3a)C A ] 



48 



+ [(3a - 449)C] + 2807>jV>Q + 192T F N f C A ] — 

+ [(1512a- 15552£(3)+ 89008) r f JV/Cj 

- 3456T F N f Cj - 84487jA^C> - l2352T$NfC A 

+ (20736^(3) + l9920)T F N f C A C F 



- (81a 3 - 162C(3)a 2 4 

- 486^(3) +75607)c]] 



162a 2 
a 3 
6912 " 



-648C(3)a + 918a 



0(a A ) 



(9) 



where a = g 2 /(16% 2 ), C,(n) is the Riemann zeta function and 
the group Casimirs are defined by 



y = Z T y (6) T a T a 



c F i , r d f bcd 



C A 8 ab , Tr(T a T b ) = T F 8 ab 
(10) 

where T" are the generators of the Lie algebra of the colour 
group. Throughout we use dimensional regularization in 
d = 4 — 2e dimensions. To establish Jy(a) we have used the 
version of the Gribov-Zwanziger Lagrangian where the y de- 
pendent terms are regarded as being part of the interaction 
Lagrangian. Therefore, the gluon and <|>?* fields are regarded 
as massless. This is the form required for applying the MIN- 
CER algorithm which determines massless 2-point functions 
to three loops and the finite part in e. Thus Jy(a) is deter- 
mined by renormalizing the Green's function where the op- 
erator associated with y 2 in (|T) is inserted at zero momentum 
into the A"-^ c Green's function. This involved 2 one loop, 43 
two loop and 1082 three loop Feynman diagrams which were 
generated with the Qgraf package, Il28ll . 

Clearly the first Slavnov-Taylor identity, ©, appears to 
hold for all a and is central to the observation that like the 
Landau gauge the other anomalous dimensions (and the p- 
function) are equivalent to their a dependent values with or 
without the localising features. By contrast the Gribov mass 
parameter is not only a dependent but is not proportional to 
the sum of the gluon and Faddeev-Popov ghost anomalous di- 
mensions as was the case in the Landau gauge. Therefore, for 
non-zero a an extra independent renormalization constant is 
required to render (Q]i finite unlike the Landau gauge. 
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III. GAP EQUATION. 



IV. INFRARED PROPERTIES. 



Having established the Gribov-Zwanziger Lagrangian as a 
tool to perform calculations one can begin to address the prob- 
lem of what its properties are at zero momentum, returning to 
the Landau gauge. First, one can construct the gap equation 
satisfied by y at two loops by evaluating the vacuum expecta- 
tion value (01. This requires the Feynman rules of ([TJ as well 
as the propagators. Due to the mixed 2-point term in (Q]) the 



A^ and ^ sector propagators are 



(AZ(p)A b v (-p)) 

(AZ(p)$ c (-p)) 



8 a V 



l(p 2 ) 2 +C A f] 



V2[(p 2 ) 2 + C A f] 



pMp) 



2 — T W 

P 

P 2 [(p 2 ) 2 + C A f] 



pAp) (ii) 



where Ppv(p) — T|,uv — PpPv/p 2 - Therefore using these to 
evaluate the 1 one loop and 17 two loop Feynman graphs con- 
tributing to (01, we arrive at the two loop MS gap equation sat- 
isfied by the running Gribov mass y, for Nf massless quarks, 
id, 



Ca 



+ 



C 



3. (C A f 
2017 11097 



768 



. 95 65 (C A r 

2048 *+ 256 « 2 ) " 48 ta ^ 



128 V \ S 



2 



1137 
2560 



V5C(2) 



205% 2 
512 



24 



a 2 + 0(a 3 ) 



(12) 



where s 2 = (2v^/9)Cl 2 (23t/3) with Cl 2 (x) the Clausen func- 
tion and p is the scale introduced in dimensional regulariza- 
tion which ensures a dimensionless coupling constant in d- 
dimensions. This calculation also required the use of the 
QGRAF package, l28Tl . to generate the electronic version of 
the Feynman diagrams as well as the use of Form to handle 
the algebra. Additionally the two loop master vacuum bubble 
integral 



I(mlmW 3 ) = l kl [k 2„ mW ^ mk _ l)2 _ ml] d3) 



where mf £ {0,/v / CaY 2 ,— iVCXj 2 } was required to the finite 
part in e. 



As Gribov was able to establish the one loop Faddeev- 
Popov ghost enhancement by evaluating the ghost 2-point 
function and applying the one loop gap equation, 11711 . it is 
now viable to extend that calculation. At two loops there are 
31 Feynman diagrams to evaluate. If we write the Faddeev- 
Popov form factor as 



D C (P 2 



1 



[l+u(p 2 )} 



(14) 



where D c (p 2 )/ p 2 represents the propagator itself, then u(p 2 ) 
represents the radiative corrections. So ghost enhancement 
then translates into the Kugo-Ojima confinement criterion of 
u(Q) = — 1, 02911 . Hence, to consider this in the Gribov- 
Zwanziger context one needs to evaluate the two loop ra- 
diative corrections in the small p 2 expansion. In practical 
terms this means carrying out the vacuum bubble expansion 
of the Feynman diagrams and using the master integral dT3b . 
for instance, and the cases where the propagators have inte- 
ger powers larger than unity. Completing this one finds that 
m(0) = — 1 at two loops precisely when y satisfies the gap 
equation of ( fT2l) . 11 111 . In extending the one loop result of 
Gribov to two loops means that the Kugo-Ojima confinement 
criterion is consistent at two loops in ([]]) as well as being a 
check on the actual explicit computation of (fT2l . Moreover, 
the fact that the condition is precisely satisfied at two loops is 
also consistent with the all orders proof given in 1 1 5 | . There 
general arguments established that the ghost propagator be- 
haves as 1 /(p 2 ) 2 as p 2 — > in (Q]). 

Repeating the analysis for the gluon propagator is more in- 
volved due to the mixing of A^ with ^ b in the quadratic part of 
the Lagrangian. If one formally writes the one loop correction 
to the 2x2 matrix of 2-point functions as, 11211 . 



p 2 8 ac 

_y2 jcab 



y2 ^acd 

-p 2 h aC h M 



( X8" c 
\ Mf cab 



j j sued \ 

N ibcd J a + 0{a 2 ) 
(15) 

where the common factor of P M v(p) in the Landau gauge has 
been factored off, N abcd = Qh ac h bd + wf ,ce f bde + Rf« be f de 



Sdt cd and 



•abed 



Tr r, 



(a b 



d) 



(16) 



where 7^' is the colour group generator in the adjoint represen- 
tation. The 14 one loop Feynman graphs contributing to the 
2-point function matrix have been evaluated both exactly and 
in the vacuum bubble expansion and are in agreement when 
the exact expression is expanded in powers of p 2 . Inverting 
the 2-point matrix gives the propagator matrix and formally 
we find the propagators to one loop are given by, lfl2Tl . 



T 8 cp 



T fPcd _ J_ftp§dq 

■> p- 



-f 



■cpq 



' {{p'Y+CAf] 

f Ab c P CfPi 



T f cdr fP1 r 



0(a 2 ) 



(17) 
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where 



V. EFFECTIVE COUPLING CONSTANT. 



q>cdpq = Gb cp b dq +Jf cpe f lcie +Kf de f pqe +Ld c A pv (18) 



and 



A = 



(-1) 



+ C A y 4 (Q + C A R + ^C A W)] . 



(19) 



The first entry in this matrix corresponds to the correction to 
the gluon propagator and in terms of the 2-point functions is 



(-1) 



+ C A y 4 (Q + C A R + \C A W)] . 



(20) 



Hence, if one writes the form factor of the gluon as D A (p 2 ) 
where D A {p 2 )P/jv(p) / p 2 is the propagator, then one finds that 
D A (p 2 ) — > as p 2 — > at one loop, 111 211 . Hence, one has 
gluon suppression at this order in (Q~|i consistent with DSE and 
lattice studies. 

We close this section on the propagators by noting that for 
the other fields one and two loop calculations reveal parallel 
behaviour. For the 03^ b field it too has a completely similar 
enhancement at two loops if one computes the two loop 2- 
point function, in a similar way to that of the Faddeev-Popov 
ghost, and applies the two loop gap equation. So the Kugo- 
Ojima confinement criterion is satisfied by both sets of anti- 
commuting fields. However, it is worth noting that the origi- 
nal Kugo-Ojima analysis was for a non-abelian gauge theory 
with only Faddeev-Popov ghosts and not with the additional 
Zwanziger ghosts. Therefore, for completeness it would seem 
that the study of [29] needs to be extended to the case of (Q]). 
Equally at one loop the exact expressions for both the co", ft and 
Faddeev-Popov ghost form factors are the same. Explicitly 
we have 



D c (p 2 ) = D a { P 2 ) 



Kp 



I in 



, , iff 



tan 



C A f 



Ap 2 
C A f 



+ 



Kp 2 ) 2 



WCaT 



- tan 



In 



P 2 

(P 2 ) 2 



c A f 



c~ A f 



C A a 



(21) 



Moreover, these facts suggest that there is an all orders equal- 
ity which, if so, could possibly be established by general sym- 
metry arguments. Also, whilst the ty^ 2-point function is 
more complicated due to the appearance of the colour group 
tensors, the channel corresponding to the tree term in the La- 
grangian is equivalent at one loop in MS to the gap equation. 
This suggests that in this channel there is a change in its in- 
frared behaviour. 



One quantity which receives wide attention in both DSE 
and lattice studies is the renormalization group invariant def- 
inition of the effective strong coupling constant based on the 
renormalization properties of the ghost gluon vertex. (See, for 
instance, H.) Since this vertex does not undergo renormaliza- 
tion in the Landau gauge, then the effective coupling constant 
is defined from the gluon and ghost form factors as 



af( P 2 ) = u(v)D A ( P 2 )(D c (p 2 )y 



(22) 



where ct(p) = g 2 (/J) / (4ji) is the running strong coupling con- 
stant. Although both numerical approaches use the same 
definition they do not definitively agree on the infrared be- 
haviour aside from the fact that there is freezing at a finite 
value. However, some DSE studies find af ff (0) = 2.97 for 
SU (3), lH^HH], whilst various lattice and DSE analyses sug- 
gest af ff (0) = 0, JUIH]. Given that we have analysed both 
form factors in the Gribov-Zwanziger Lagrangian, it is a sim- 
ple exercise to evaluate d22b in the p 2 — > limit. We find that 
at one loop in the MS scheme, 11211 . 



ocf ff (0) = lim 



aOi)[l+C A (|L-||)a] (p 



C A f[l+C A ^ 



L — 



where L = In (^) . Thus for St/ (3) we find af ff (0) 



(23) 
1.768 



independent of the number of massless quarks. Whilst this is 
less than the value for the DSE approach it is finite and non- 
zero. Moreover, the N c dependence is in qualitative agreement 
with 03(X 13111 . In the determination of a finite value it is worth 
noting the key role played by the gap equation in ensuring the 
cancellation of the dimensionful and dimensionless quantities 
to leave a pure number. We also remark that the C0° 6 gluon 
vertex has a similar freezing to the same value, fl2il . 

One property of ( 1221 which motivated recent studies of Lan- 
dau gauge QCD and the suggested condensation of low di- 
mensional operators was the numerical evidence of a 1 lp 2 
power correction to the effective coupling constant, l34l l35i 



]. Although that analysis was in the MOM scheme, if 
it is a more general property of the Landau gauge effective 
coupling constant it ought to be present in other schemes. As 
we have exact expressions for the form factors it is possible 
to evaluate ( 1221 in powers of p 2 and see whether a l/p 2 term 
emerges as a correction or if the expected l/(p 2 ) 2 term is the 
first term. The latter is certainly expected to be present due to 
the existence of the gauge invariant dimension four operator 
(G^ v ) 2 which can condense. Therefore, returning to ( TT3T > and 
formally expanding in powers of y 2 , we have at leading order 
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at one loop, I12H . 



97 \ 2 3tcvW 
36j^ + 8 
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T F Nfp 2 
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M 



U 



IIC A 



r 



Q = 

W = R 



4 \p 

,2 



C A 



S = 0{a l ) . 
Hence, at one loop 

3C A K y/Ckf 



Da( P 2 ) 
D C (P 2 ) 



-1 



4 \v 2 



Caci 



(24) 



(25) 



-l 



Thus evaluating ocf (j> 2 ) in the same expansion and identify- 
ing the y independent pieces as deriving from the perturbative 
part of the coupling constant, we find at one loop 



«t f V 



afV) 



9C A n VClf 



O 



(P 



2\2 



(26) 

Thus a 1/j? 2 correction emerges and the power correction of 
i34l l35l 13a, l37ll can be mimicked by the presence of the 
Gribov mass in the Gribov-Zwanziger Lagrangian. Though 
this ought to be qualified by noting that the condensation of 
the gauge variant operator ^(A") 2 has also received a large 
amount of attention in recent years in the Landau gauge, 
lHHHIl], and could be an alternative justification for the lattice 
result. It is also worth observing that both y 2 and 5 (A") 2 have 
the same anomalous dimensions in the Landau gauge M40ll4lll . 
Indeed the effect that the Gribov mass has on the condensation 
of j(A") 2 has been examined in 11221 l42l I43TI . One final point 
is that the sign of y 2 is not fixed in (Q~|i. Therefore, one can alter 
the sign of y 2 in djot to obtain the same sign as that in ll36[[3~7ll . 
Although this is for a different scheme (MOM) compared to 
our MS result, it would result in an effective tachyonic gluon 



mass. This may appear unacceptable. However, as the gluon 
is not truly an observable field due to confinement this ought 
not to be a hindrance. More specifically this tachyonic nature 
is only revealed in the next to high energy limit. If one ex- 
amines the full gluon propagator at one loop, there is no pole 
in the propagator. Therefore, one cannot truly identify a mass 
for the confined gluon in the conventional way that one does 
for a fundamental unconfined field. Another observation from 
a phenomenological point of view is that in ll44ll it was noted 
that a tachyonic gluon mass in current correlators could repro- 
duce experimental data more accurately than either a massless 
gluon or one with a non-tachyonic mass. 

VL CONCLUSION. 

We have studied the Gribov-Zwanziger Lagrangian, flill . 
at one loop and established that it qualitatively reproduces 
the propagator and effective coupling constant behaviour ob- 
served in the infrared using DSE and lattice studies. More- 
over, the gap equation has been determined at two loops in 
MS and ensures the Kugo-Ojima confinement criterion is sat- 
isfied. This would appear to establish the Gribov-Zwanziger 
Lagrangian as a useful tool to examine other phenomenolog- 
ical aspects of QCD such as the operator product expansion 
and the evaluation of vacuum expectation values of gauge in- 
variant operators, such as (G^ v ) 2 . Further, we have examined 
the renormalization of the Gribov-Zwanziger Lagrangian in 
linear covariant gauges albeit in the limit of small gauge fixing 
parameter. Clearly it would be a formidable task to fully ex- 
tend the Landau gauge Gribov-Zwanziger Lagrangian to other 
gau ges such as linear covariant and maximal abelian gauges, 
B241 145I1 . However, if it were possible then it would be in- 
teresting to extend our present Landau gauge propagator and 
effective coupling constant analysis to these gauges in order to 
understand what their infrared structure is. This would com- 
plement DSE and lattice analyses. 
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